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ABSTRACT 

The location of M dwarfs in the V — K s -Mks color-magnitude diagram (CMD) has been shown to correlate with 
metallicity. We demonstrate that previous empirical photometric calibrations of M dwarf metallicity exploiting this 
correlation systematically underestimate or overestimate metallicity at the extremes of their range. We improve upon 
previous calibrations in three ways. First, we use both a volume-limited and kinematically-matched sample of F and G 
dwarfs from the Geneva-Copehnagen Survey (GCS) to infer the mean metallicity of M dwarfs in the Solar Neighborhood. 
Second, we use theoretical models of M dwarf interiors and atmospheres to determine the effect of metallicity on M 
dwarfs in the V — K s -Mks CMD. Third, though we use the GCS to infer the mean metallicity of M dwarfs in the 
Solar Neighborhood, our final calibration is based purely on high-resolution spectroscopy of FGK primaries with M 
dwarf companions as well as the trigonometric parallaxes and apparent V- and its-band magnitudes of those M dwarf 
companions. As a result, our photometric calibration explains an order of magnitude more of the variance in the 
calibration sample than previous photometric calibrations. We use our calibration to non-parametrically quantify the 
significance of the observation that M dwarfs that host exoplanets are preferentially in a region of the V — K s -Mks 
plane populated by metal-rich M dwarfs. We find that the probability p that planet-hosting M dwarfs are distributed 
across the V — K s -Mks CMD in the same way as field M dwarfs is p = 0.06 ± 0.008. Interestingly, the subsample of M 
dwarfs that host Neptune and sub-Neptune mass planets may also be preferentially located in the region of the V — K s - 
Mks plane populated by high-metallicity M dwarfs. The probability of this occurrence by chance is p = 0.40 ± 0.02, 
and this observation hints that low-mass planets may be more likely to be found around metal-rich M dwarfs. The 
confirmation of this hint would be in contrast to the result obtained for FGK stars, where it appears that metal-rich 
and metal-poor stars hosts Neptune-mass planets with approximately equal probability. An increased rate of low-mass 
planet occurrence around metal-rich M dwarfs would be a natural consequence of the core-accretion model of planet 
formation. 

Key words, planetary systems — planets and satellites: formation — stars: abundances — stars: low-mass — stars: 
statistics 



1. Introduction 

The determination of metallicity for M dwarfs is a very 
difficult problem (e.g. lGustafssonl ll989). Their cool atmo- 
spheres permit the existence of many molecules for which 
molecular opacities are currently poorly constrained. As a 
result, the estimation of the continuum level of a spectrum 
is challenging, rendering line-based metallicity indicators 
unreliable. The poorly constrained molecular opacity data 
currently available makes the determination of metallicity 
through spectral synthesis difficult as well. For those rea- 
sons, alternative methods must be employed to estimate M 
dwarf metallicities. 

The main sequence lifetimes of M dwarfs are longer than 
the Hubble time, so they have not yet departed much from 
the zero-age main sequence. Consequently, M dwarfs might 
be expected to form a two-parameter sequence in mass and 
metallicity suggesting that a two-color broad-band photo- 
metric calibration might constrain their properties. There 
have been several attempts to obtain the metallicity of M 
dwarfs using their p hotometric properties , i nclud ing two re- 
cent breakthroughs. iBonfils et al.l (|2005af) — IB05I hereafter - 
had the subtle insight to realize that M dwarfs in binary or 



multiple systems should have metallicities commensurate 
with the easily-measure d me tallicity of an FGK primary in 
the system. In that way. IB05I identified a calibration sample 
of M dwarfs with metallicities securely determined in one of 
two ways: (1) high-resolution spectroscopy of an FGK com- 
panion and (2) high-resolution spectroscopy of M dwarfs for 
which T e f f and log g could be fixed with photometric data. 
For the former, the metallicity is very likely the same as 
the metallicity inferred from high-resolution spectroscopy 
of its FGK companion. For the latter, spectral synthesis af- 
ter fixing T e ff and \ogg with photometric data eliminates 
some degeneracy and produces a reasonable metallicity es- 
timate. They noted that low-metallicity M dwarfs have blue 
V— K s color at constant if s -band absolute magnitude Mk s , 
and they fit a linear model to their calib ration sample us- 
i ng V — K s an d Mk s to predict [Fe/H]. I Johnson "fcA pps 
(2009) - JA09 hereafter - addressed the relative la ck of 
high-metallicity M dwarfs in the calibration sample of IB05I 
and created an empirical model in which the distance of 
an M dwarf above the field M dwarf main sequence (MS) 
in the V — K s -Mk s color-magnitude diagr am (C MD) in- 
dicated its metallicity. The great insight of JA09 was that 
the mean metallicity of a population of M dwarfs could be 
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characterized by the easily-measured mean metallicity of 
a similar population of FGK stars. Indeed, they assumed 
that the field M dwarf MS was an isometallicity contour 
with the same metallicity as a volume-limited sample of G 
and K stars and fit a linear model using the distance above 
the field M dwarf MS to predict [Fe/H]. 

M dwarfs are attractive targets around which to search 
for low-mass planets because they have large reflex ve- 
locities and transit depths even for low-mass and small- 
radius companions. Given that the metallicity of proto- 
planetary disks is a key parameter in models of p lanet 
formation (e.g. lLaughlin et al.ll200i llda fe Lhl 12004) . the 
metallicity of M dwarfs that host planets will constrain 
the planet formation process in low-mass protoplane- 
tary disks. Indeed, it is well-established that metal-rich 
FGK stars are more likely to host giant planets (e.g. 
ISantos et ahl 12004 iFischer fe Valentil 120051 ). but there is 
also evidence to suggest that metal-rich FGK stars are 
not much more likely to host Neptune- mass planets than 
their low-metallic ity counterparts (e .g. lUdrv et al.l [2 006: 
ISousa et al]|2008HBouchv et al.ll2009[ ). Already. UA09I have 
used their model of M dwarf metallicity to suggest that the 
M dwarfs that host planets are preferentially metal-rich. 
However, UA09I did not address whether the apparent lack 
of a correlation between FGK host stellar metallicity and 
the presence of Neptune-mass planets extends to M dwarfs. 

In this paper, we examine a calibration sample of M 
dwarfs with secure ly es tima ted m etallicities and we show 
that the models of B05l and UA09I systematically underes- 
timate or overestimate metallicity at the extremes of the 
range of this calibration sample. We demonstrate that a 
volume-limited and kinematically-matched sample of Sun- 
like stars produces a better estimate of the mean M dwarf 
metallicity in the Solar Neighborhood, and we use M dwar f 
models of different metallicities from iBaraffe et al.l (jl998) 
to improve on the technique described in The posi- 

tion of an M dwarf in the V — K s -Mk s CMD remains an 
indicator of its metallicity, and we use that fact to non- 
parametrically quantify the significance of the observation 
that planet-hosting M dwarfs are preferentially in a region 
of the V— K s -Mk s plane populated by metal-rich M dwarfs. 
Moreover, we identify for the first time a hint that the sub- 
sample of M dwarfs that host Neptune and sub-Neptune 
mass planets may also be more likely to be in the region 
of the V — K s -Mks CMD associated with metal-rich M 
dwarfs. We describe our analysis in §2 and summarize our 
findings in §3. 



2. Analysis 



2.1. Testing Previous Calibrations 



We first collect from lB05l and lJA09l a calibration sample of 
M dwarfs in wide binary or multiple systems with an FGK 
primary. The metallicity of the FGK primary in the system 
is straightforward to measure from a high-resolution spec- 
trum, and if the M dwarf secondary and the FGK primary 
formed in the same molecular core, then the expectation is 
that the two should ha ve co mmensurate metallicities. We 
collect 13 examples from lB05l selecting only those M dwarfs 
with precise T^-band magnitudes from CCD photo metry . 
We also collect six high-metallicity examples from lJA09i . 
We summarize this calibration sample of M dwarfs in the 
first 19 lines of Table [U 



We compute the met allicit y p redicte d for this calibra- 
tion sample from both the lB05l andlJA09 relati ons a nd com- 
pare it to the observed values. Note that the IB05I relation 
was initially based on a calibration sample that included M 
dwarfs in binary or multiple systems with FGK primaries 
(some with V^-band magnitudes deduced from photographic 
plates) as well as low-metallicity M dwarfs with metallicity 
inferred from spectroscopy after fixingJT e / f and log g with 
photometric data. Meanwhile, the |JA09 calibration sample 
included only the six metal-rich M dwarfs in binary or mul- 
tiple systems with FGK primaries listed in rows 14 through 
19 in Table [TJ For those reasons, we believe that the 19 
M dwarfs with metallicities inferred from high-resolution 
spectroscopy of FGK primaries in Table [T] is the largest 
and most reliable set of M dwarf metallicities fro m w hich 
to ver ify previous calibrations. We apply both the lB05l and 
IJA09I relations to this sample and compute the residual 
between each model and observation. We plot the distri- 
bution of re sidua ls for both models in Figure [TJ and we 
find that the lB05l relation syst ematic ally underestimates M 
dwarf metallicity and that t hel.lAO!) relation systematically 
overestimates M dwarf metallicity. 

2.2. A Physically-Motivated Empirical Model of M Dwarf 
Metallicity 



As di scussed in §2.1 and Figure [H the models of [BOl and 
JA09I have non-negligible residuals when applied to the cal- 
ibration sample in Table [T] Still, there is a correlation be- 
tween the metallicity of an M dwarf and its distance in 
the V - K s -M Ks plane f rom t he field M dwarf MS. We 
attempted to improve the JA09| model by reassessing both 
the zero point of the model and the direction from the M 
dwarf MS in the V — K s -Mk s plane best correlated with 
metallicity. 

Recall that IJA09I set the mean metallicity of the Solar 
Neighborhood M dwarf sample equal to the mean metal- 
licity of a volume-limited sample of G 0-K2 stars (4.0 < 
My < 6.5) from the SPOCS catalo g of IValenti fe Fischerl 
(|2005l ). The SPOCS sample of lJA09l was based on a catalog 
of stars selected to have absorption lines deep enough to en- 
able high-precision radial velocity detection of exoplanets. 
As a result, the SPOCS sample is biased against metal-poor 
stars and therefore potentially unsuitable for the determi- 
nation of the average metallicity in the Solar Neighborhood 
(as noted in IValenti fe Fischerl 12005). Alternatively, the 



(a,s notea m ivaienti & lus cncr zuuo). Alternatively, tne 
Geneva-Copenhagen Survey CGCS - iNordstrom et alJl200l 
IHolmberg et~alll2007l 120091 ) of Solar Neighborhood F and 
G dwarfs is magnitude-complete, kinematically-unbiased, 
and free of the line depth bias inherent in the SPOCS cat- 
alog. Though the GCS metallicity estimates are based on 
Stromgren uvby/3 photometry and not high-resolution spec- 
troscopy, the precision of the GCS metallicities are suffi- 
cient when combined with the reduced bias of the sample to 
provide a better estimate of the mean Solar Neighborhood 
metallicity than the SPOCS sample. 

In addition, the UVW kinematics of a volume-limited 
sample of M dwarfs is not necessarily equivalent to the kine- 
matics of a volume-limited sample of FGK dwarfs. Since the 
mean of a sample is sensitive to outliers, and because Sun- 
like stars with outlier kinematics are also likely to be out- 
liers in metallicity, a kinematic-match is important to deter- 
mine the mean metallicity of Solar Neighborhood M dwarfs 
in this way. To address this point, we use the M dwarf UVW 
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distribution described by lHawlev et al.l (|l996f) to create a 
volume-limited and kinematicall y -matc hed sample of F and 
G dwarfs from iHolmberg et al.l ([20091 ) from which we in- 
fer the average metallicity of the Solar Neighborhood M 
dwarf population. In Figure [21 we superimpose the UVW 
velocity-space distri bution of local M dwarfs derived by 
Haw lev et al.l (|l996h on top of the UVW velocity-space 
distribution of F and G stars from the GCS with parallax- 
based distance estimates that place them within 20 pc of 
the Sun. We bootstrap resample from the subset of GCS 
stars within 20 pc and with kinematics c onsistent with 
the M dwarf velocity ellipsoid as defined in lHawlev et al.l 
(|!996t) . We ensure that 68% of the GCS stars in each 
bootstrap sample have UVW v elocities that place them 
within the one-sigma contour of lHawlev et al.l (jT996) and 
that the rest of each bootstrap sample lies within the two- 
sigma contour. In the end, we find that a volume-limited 
and kinematically-matched sample of F and G dwarfs 
from the GCS survey has a mean metallicity of [Fe/H] 
= —0.14 ±0.06. We obtain a similar result with a sample of 
GCS stars volume-limited in the same way as the volume- 
limited SPOCS sample of lJA09l . for which we find a mean 
Solar Neighborhood metallicity of [Fe/H] = -0.15 ± 0.02. 
In this case, the superior statistics of the larger volume- 
limited sample is enough to formally achieve a higher pre- 
cision than the volume-limited and kinematically-matched 
sample, though the volume-limited sample is subject to a 
greater degree of possible systematic error. For that reason, 
we regard the volume-limited and kinematically-matched 
result as likely more reliable. 

The mean metallicity of our volume-limited and 
kinematically-matched GCS sa mple suggests that the field 
M dwarf MS defined by UA09I is an isometallicity contour 
with [Fe/H] w —0.14. Note that if the iso metal licity con- 
tour corresponded to [Fe/H] — —0.05 as in lJA09l . five stars 
from Table [1] with [Fe/H] < -0.05 would be to the right 
of the isometallicity contour indicating [Fe/H] > —0.05. 
Alternatively, if we assume that the isometallicity contour 
corresponds to [Fe/H] w —0.14, then only two of the 19 
stars are on the wrong side of contour. 

We now determine which direction in the V — K s -Mk s 
plane an isochrone moves as a function of metallicity. In 
Figure [3] we plot the M dwarfs with securely determined 
meta llicity from Table [TJ along with the M dwarf MS 
from | JA09l and two different isochrones from IBaraffe et all 
(|1998h . We use the transforma tion of Carpenter! (l200il ) 
to transform the Kqit given in IBaraffe et al.l (|l998h into 
K s . The left- most isochrone corresponds to a population 
with [Fe/H] = -0.5 and Y = 0.25 while the right-most 
isochrone corresponds to a population with [Fe/H] = and 
Y = 0.275. Both isochrones use mixing-length parameter 
I = 1 for a 5 Gyr population (there is no detectable evo- 
lution in V - K s -M Ks CMD after 3 Gyr). The horizontal 
lines connect points of constant mass. With all other pa- 
rameters constant, metallicity should best correlate with 
horizontal shifts in the V — K s -Mk s plane. For that rea- 
son, wc compute the distance from the M dwarf MS in the 
horizontal direction for each M dwarf with secure metal- 
licity from Table [TJ We then fit a linear model using this 
distance as a predictor with [Fe/H] as the response. We find 
that 



A(V- K s ) ee (V - K s ) obs -(V- K s ) lso 

is the optimal model. In this case, Mk s as func- 
tion of V — K s is given by the fifth-order polynomial 
with coefficients in increasing order (- 9.5893 3,17.3952,- 
8.88365,2.22598,0.258854,0.0113399) from[TA0l To aid in 
the calculation of (V — K s ) iso , we give the same curve with 
V — K s as a function Mr- s : it is a fifth-order polynomial with 
coefficients in increasing order (51.1413,-39.3756,12.2862,- 
1.83916,0.134266,-0.00382023). 

We use tw o of the model selection criteria given in 
iHockind (|1976f) to evaluate all three models. First, we com- 
pute the residual mean square (RMS), defined as 



RMS r 



SSE 



n — p 



(2) 



where n is the number of data points, p is the number 
of predictors in the model, and SSE p is the residual sum 
of squares for a p-term model. In general, models with 
smaller values of RMS P are best-suited to prediction. For 
our m odel, we find that RMS P = 0.02; for the mod el of 
UA09I the value is RMS P = 0.04 while for the model of lB05l 
the value is RMS P = 0.05. Next, we compute the adjusted 
square of the multiple correlation coefficient Ra P , which is 
widely used to judge the fit of a model. A value of R 2 p = 1 
indicates that a model explains all of the variance in a sam- 
ple, while R\ p = indicates that the model explains none 
of the variance. R 2 ap is defined as 



R 2 ap = l-i/,- i) 
SST 



RMS r 



" K " SST 



(3) 



[Fe/H] = 0.79A (V - K s ) - 0.17 



(1) 



where yi and y are the sample and its mean, respectively. 
For our model, we find that R? ap = 0.49; for the model of 
UA09I the value is R 2 ap = 0.059 while for the model of lB05l 
the value is R 2 p < 0.05. We note that our model explains 
almost an order of magnitude more of the varianc e in t he 
calibr ation sample than either model presented in IB05I or 
I.TA09L 

Differing T e ff scales have been well-noted as a source 
of metallicity discrepan cies in metallicity st udies of the 
Solar Neighborhood fe.g. lHolmberg et al.ll2007t ). The differ- 
ing T e ff calibrations between the GCS and other surveys 
will not affect our results, as our calibration (including the 
metallicity of the M dwarf MS) is based on the horizon- 
tal distance i n the V — K s -Mk s CMD from the mean M 
dwarf MS of I JA09I of M dwarfs with metallicities known 
from high-resolution spectroscopy of FGK primaries. We 
only used the GCS Stromgren-based metallicities to estab- 
lish the fact that the mean Solar Neighborhood metallicity 
is closer to [Fe/H ] « -0.15 than it is to [Fe/H] « -0.05 as 
arg ued by[jA09. For that reason, our calibration is based 
purely on high-resolution spectroscopy of FGK primaries as 
well as the trigonometric parallaxes and apparent V- and 
ifs-band magnitudes of their M dwarf companions. 

The M dwarfs in binary systems with FGK primaries 
that we use to fix our calibration are not a volume-limited 
or kinematically-matched sample. The volume-limit and 
kinematic-match were only necessary to verify the fact that 
the mean metallicity of the M dwarf population in the Solar 
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Neighborhood is a well-defined quantity. That verification 
is a nece ssary precondition that must be established be- 
fore any I JA09I style calibration using distance from the 
field M dwarf MS in the V - K s -M Ks CMD can even be 
considered. Once the points along the field M dwarf MS 
in the V — K s -Mk s CMD are fixed to the mean metal- 
licity of the Solar Neighborhood, the metallicity of an M 
dwarf with given V — K s color and absolute magnitude 
Mks along the curve is specified regardless of its position 
or velocity. Indeed, when we build our calibration using 
only the metallicities of M dwarfs in binaries with FGK pri- 
maries, their V- and K s -band magnitudes, and trigonomet- 
ric parallaxes, we find that the mean metallicity of the Solar 
Neighborhood M dwarf MS based on the calibration sam- 
ple ([Fe/H] = — 0.17±0.07) is statistically indistinguishable 
from the mean metallicity inferred from the volume-limited 
and kinematically- matched sample ([Fc/H] = — O.I4±0.06). 

2.3. The Metallicity of M Dwarfs that Host Planets 

We plot the location of M dwarfs that host planets in 
the left-hand panel of Figure 2] The M dwarf models of 
iBaraffe et al.l (|1998l ) suggest that horizontal distance in the 
V — K s -Mks CMD best correlates with metallicity. To 
non-parametrically determine the degree to which planet- 
hosting M dwarfs are preferentially found to the right of 
the M dwarf MS, we need to quantify the likelihood that 
the cumulative horizontal distance from the isometallicity 
contour of a randomly selected sample of field M dwarfs 
can be as large as that observed in the sample of M dwarfs 
that host planets simply by chance. 

To address this issue, we create a cont rol sample of 
field M dwarfs selected from the Hipparcos (Ivan Leeuwenl 
12007ft and Yale Parallax Catalog s (jvan Altena et al.lll995l) . 
We include in the control sample those M dwarfs from the 
Hipparcos catalog that have parallaxes n > 100 mas precise 
to better than 5% and those M dwarfs from the Yale catalog 
that have parallaxes it > 100 mas. We use F-band photom- 
etry from each catalog and we obtain K s photometry for 
both samples from the 2MASS database (jSkrutskie et al.1 
l2006h . For a sample of size n we can compute the statistic 
E: 



n 

E = £ (V - K s \ -(V- K s ) iso (4) 

i=l 

To characterize the likelihood that an observed value of E 
for a subsample with size m of M dwarfs that host planets 
can be produced by chance, we use a Monte Carlo simula- 
tion. We randomly select a sample of size m from the 127 M 
dwarfs in the control sample, compute E for that sample, 
save the result, and repeat the calculation 1000 times. In 
that way, we can determine the distribution of E expected 
under the null hypothesis that M dwarfs that host plan- 
ets are distributed in the V — K s -Mk s CMD in the same 
way as field M dwarfs. We consider three sub-samples: (1) 
all planets hosts, (2) hosts of Jupiter-mass planets, and (3) 
hosts of Neptune-mass (and below) planets. We find that in 
case (1) E = 3.43 indicating only ap = 0.06 ± 0.008 proba- 
bility that the cumulative distance of the sample from the 
isometallicity contour occurred by chance. In case (2) we 
find E = 2.39 => p = 0.04 ± 0.005 and in case (3) we find 



E = 1.04 => p = 0.40±0.02. We summarize this calculation 
in the right-hand panel of Figure [4] 

2.4. Discussion 

The apparent position of planet-hosting M dwarfs in the 
region of the V — K s -Mk s CMD associated with known 
high-metallicity M dwarfs tentatively suggests that metal- 
rich M dwarfs are more likely to host Jupiter-mass and 
possibly Neptune-mass planets as well. If this correlation is 
confirmed in the future, it can be understood as a natural 
cons equence of the core- accretion model of planet formation 
(e.g. lLaughlin et al.ll200l llda fc Linll20"ol . Indeed, a more 
metal-rich protoplanetary disk will almost certainly have a 
higher surface density of solids, and that increased surface 
density enables the rapid formation of the several Earth- 
mass cores necessary to accrete gas from the protoplane- 
tary disk before the gaseous disk is dissipated. Moreover, it 
would be especially interesting if the tentatively suggested 
correlation extends to the hosts of Neptune-mass planets, as 
current evidence seems to suggest that the probability that 
an FGK star hosts a Neptu ne-mass planet is not a strong 
function of metallicity (e .g. lUdrv et all 120061 : ISousa et al.l 
[2001 IBouchv et al.ll2009h . 

If the tentatively suggested correlation between the 
presence of planets and the metallicity of their host M 
dwarfs is eventually confirmed, it might indicate a lower- 
limit on the amount of solid material necessary to form 
planets. To see why, recall that the mass of a protoplane- 
tary disk scales roughly as Mdisk oc M* and that the frac- 
tion of solid material in a disk f so iid scales roughly as 
where is the metal content of the host star. The total 
amount of solid material in a protoplanetary disk will then 
scale like M soHd cx J S oiidM dlsk oc Z»M*. Minimum-mass 
Solar Nebula models (MMSN - lHavashil fl98ll) and obser- 
vations of T Tauri disks in star-forming regions suggest 
that protoplanetary disks around young Solar-type stars 
are about 1% the mass of their host stars , albeit with sig- 
nificant scatter (e.g. lHartmann et al.lll998l ). Combined with 
the fact that the metal content of the Sun is Zq — 0.0176 
by mass, the total solid mass in the MMSN was about 
M solid ~ 60 M®. This is a lower-limit, as more careful 
calculations suggest that the proto planetary disk around 
the Sun had M soUd ~ 100 M ffi (e.g. lLissauenll993l) . In ei- 
ther case, this is a factor of a few to ten greater than the 
10 M s of material necessary to form the core of gas or ice 
giant planet in the core-accretion model of planet forma- 
tion. In the case of a Solar-metallicity mid-M dwarf with 
A/* = 0.3 Mq, the total amount of solid material in the 
disk is 70% less, about M so ud ~ 20 M e . This is factor of 
order unity to a few times the mass necessary to form the 
core of a gas or ice giant. Since planet formation likely does 
not lock-up the entire solid component of a protoplanetary 
disk in planets, reducing the total mass of solids in the disk 
- either by reducing the metallicity or mass - will also re- 
duce the chances of forming a 10 Mq core (and therefore 
a gas or ice giant) before the parent protoplanetary disk is 
dissipated. 

The confirmation of the hint of a correlation between 
the presence of low-mass planets and M dwarf metallicity 
could be evidence of this threshold solid mass necessary to 
form Neptune-mass planets. Around FGK stars, the same 
threshold solid mass suggests that a correlation between the 
presence of low- mass planets and host star metallicity might 
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occur at one-third Solar metallicity, or [Fc/H] = —0.5. This 
is just below the typical metallicity of stars observed at 
high radial velocity precision with HARPS ()Sousa et al.l 
2008) . This expected correlation might be verified as larger 
samples of low-metallicity stars are surveyed at high ra- 
dial velocity precision or by transit surveys of nearby low- 
metallicity open clusters (e.g. NGC 752 or IC 4756). 

The hint of a correlation between the presence of 
Neptune-mass (and below) planets and M dwarf metal- 
licity tentatively suggests that searches for low-mass 
planets around M dwarfs lik e the MEarth Project 
(|Nutzman fc Charbonnea"ull2008h could improve their yield 
by shading their target list toward M dwarfs that have 
red V — K s colors at constant i^ s -band absolute magni- 
tude Mks- Note that since the absolute magnitude Mk s 
of an M dwarf depends only on the logarithm of its often 
poorly-known distance, while V — K s depends linearly on 
its often poorly-known y-band magnitude, the collection of 
high-quality CCD-based T^-band magnitudes for M dwarfs 
in the Solar Neighborhood could be the first step towards 
maximizing the yield of planets around M dwarfs. 

3. Conclusion 

We showed that previous empirical photometric calibra- 
tions of M dwarf metallicity systematically underestimate 
or overestimate metallicity at the extremes of their range. 
We derived a physically-motivated model that explains 
an order of magnitude more of the variance i n the cal- 
ibration sample t hat e ither the iBonfils et al.l (|2005al) or 
iJohnson fe AppsI (|2009t ) models. We used the correlation 
underlying our model to non-parametrically show that the 
probability p that there is no relationship between position 
of an M dwarf in the V — K s -Mk s CMD and the presence 
or absence of planets is p = 0.06 ±0.008. For the subsample 
of M dwarfs that host Jupiter-mass planets, the probability 
that there is no correlation is p — 0.04 ± 0.005. Meanwhile, 
for the subsample of M dwarfs that host Neptune-mass 
(or below) planets, we find that the probability that there 
is no correlation is p = 0.40 ± 0.02. Since the models of 
iBaraffe et al.l (|1998h suggest that the position of an M dwarf 
in the V — K s -Mk s CMD is a qualitative indicator of metal- 
licity, this observation tentatively suggests that metal-rich 
M dwarfs are more likely to host planets and hints that the 
correlation may extend to low-mass planets as well. If this 
correlation is confirmed in the future, it will be in contrast 
to planetary systems around FGK stars, in which there ap- 
pears to be only a weak connection between metallicity and 
the presence of Neptune-mass planets. 
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F ig. 1. Optimally-smoothed residual distributions for lBonfils et al.l ()2005af) — IB05I hereafter - and I Johnson fc Appsl ([2009) 
JA09I hereafter. In both cases the vertical dashed line indicates the mean of the distribution. The mean value of the 



residuals is 0.12 with standard deviation 0.16, while the mean value of the I JA09I residuals is -0.12 with standard 



deviation 0.12. Note that the B05l distribution has a heavy tail at large positive values (indicating systematically low 
[Fe/H] estimates) and the IJA09T distribution has a heavy tail at large negative values (indicating systematically high 
[Fe/H] estimates). 



6 



Schlaufman & Laughlin: M Dwarf Metallicity 



^^^^^^^ 



^— M Dwarf Velocity Ellipsoid from Hawley et al. (1 996) 
FGK Stars from Holmberg et al. (2009) 
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Fig. 2. Velocity ellipsoids inferred for a volume-limited sample of early M dwarfs from lHawlev et all ()1996l ) superimposed 
on the WW velocit y distribution of a volu me-limited sample (d < 20 pc) of Sun-like stars from the Geneva-Copenhagen 
Survey (gray points - iHolmberg et al.ll2009h . The light curve denote s the o ne-sigma region while the heavy curve denotes 
the two-sigma region. Bootstrap resampling of the IHolmberg et"aT1 ()2009f ) sample with the constraint that 68% of each 
bootstrap sample lies within the one-sigma contour and that the other 32% lies within the two-sigma contour produces 
a volume-limited and kinematically-matched population metallicity of [Fe/H] = —0.14 ± 0.06. 
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Fig. 3. Position of M dwarfs with se cure metallicities from Table [T] (blue points) in the V — K s -Mk ? CMD in relation 
to the field M dwarf MS from JA09 (black line) and the theoretical isochrones of iBaraffe et al. (1998) (blue lines). The 
color of the isochrone line gives its metallicity: [Fe/H] = —0.5 and Y — 0.25 on the left and [Fe/H] = and Y — 0.275 
on the right. Both isochrones assume mixing-length parameter I = 1 for a 5 Gyr population, as there is no detectable 
evolution in the V — K s -Mk s CMD after after about 3 Gyr. The horizontal lines connect points of constant mass. The 
models indicate that differences in metallicity should best correlate with horizontal shifts in the V — K s -Mk s plane. 
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Fig. 4. Left: Position of M dwarfs known to host Jupiter- mass planets (dark blue triangles) and Neptune-mass (or below) 
plane ts (blue squares) in relation to a control sample of field M dwarfs (gray points) and the field M dwarf MS from 
I JA09I (black line). Again, like the high-metallicity M dwarfs, the M dwarfs that host planets are concentrated to the 
right of the field M dwarf MS. Right: Distribution of cumulative sample distances from the field M dwarf MS of JA09, 
which we assume to be a [Fc/H] « —0.17 isometallicity contour in the V — K s -Mks CMD. Points to right of the field 
M dwarf MS add their distance to the sum, while points left of MS subtract their distance from the sum. We generate 
each distribution with a Monte Carlo simulation. First, we randomly select from the field M dwarf sample a number 
of stars equal to the number of M dwarfs known to host planets of a certain type. We then compute the cumulative 
horizontal distance of that random subsample from the field M dwarf MS. We repeat this process 1000 times to generate 
the distribution of sample cumulative horizontal distances from the held M dwarf MS given no correlation between the 
presen ce of an exoplanet and the location of its host in the V — K s -Mk s CMD. In all cases, we confirm the findings of 
I JA09I that the M dwarfs that host exoplanets are preferentially to the right of the field M dwarf MS. In particular, we 
find that the probability p that there is no correlation between the location of an M dwarf in the V — K s -Mk s CMD 
and its status as an exoplanet host is p = 0.06 ± 0.008. For the subsample that hosts Jupiter-mass planets, we find that 
the probability is p = 0.04 ± 0.005. More interestingly, we find that the probability that there is no correlation between 
the location in the V — K s -Mk s CMD and an M dwarf's status as the host of a Neptune-mass (or bel ow) exoplanet 
is p = 0.40 ± 0.02. If M dwarfs to the right of the field M dwarf MS are metal-rich as suggested by the iBaraffe et al.l 
(|1998|) models and argued bv lJA09l . then this observation may be evidence for an increased incidence of low- mass planets 
around metal-rich low-mass stars, a trend which is not observed in FGK stars. 
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Table 1. M Dwarfs in Binary Systems with an FGK Primary and Those That Host Planets 



System 


V 


K s 


7T 


V -K a 


M Ka 


Spectroscopic 


Photometric 


Comment 


Reference 






[mag] 


[mas] 


[mag] 


[mag] 


[Fe/H] 


[Fe/H] 






Gl 105B 


11.67 


6.57 


138.72 


5.1 


7.28 


-0.19 


-0.14 




1,2,3,4 


Gl 107B 


10.06 


5.87 


89.03 


4.19 


5.62 


-0.03 


0.01 




1,2,4,5 


Gl 166C 


11.17 


5.9 


198.08 


5.27 


7.38 


-0.33 


-0.05 




1,2,4,5 


Gl 212 


9.8 


5.76 


80.13 


4.04 


5.28 


0.04 


0.09 




1,2,4,6 


Gl 231. IB 


13.42 


8.28 


51.76 


5.14 


6.85 


-0.02 


0.12 




1,2,4,5 


Gl 250B 


10.09 


5.72 


114.94 


4.37 


6.02 


-0.15 


-0.07 




1,2,4,5 


Gl 297.2B 


11.8 


7.42 


44.47 


4.38 


5.66 


-0.09 


0.13 




1,2,3,4 


Gl 324B 


13.14 


7.67 


76.8 


5.47 


7.10 


0.32 


0.26 




1,2,4,5 


Gl 53. IB 


13.6 


8.67 


48.2 


4.93 


7.09 


0.07 


-0.17 




1,2,4,5 


Gl 768. IB 


13.1 


8.01 


51.57 


5.09 


6.57 


0.07 


0.22 




1,2,4,5 


Gl 783.2B 


13.94 


8.88 


48.83 


5.06 


7.32 


-0.16 


-0.19 




1,2,4,5 


Gl 797B 


11.88 


7.42 


47.65 


4.46 


5.81 


-0.07 


0.11 




1,2,4,5 


Gl 81. IB 


11.21 


7.6 


29.43 


3.61 


4.94 


0.09 


-0.08 




1,2,4,5 


HD 46375B 


11.8 


7.84 


29.93 


3.96 


5.13 


0.24 


0.12 




2,4,7 


HD 38529B 


13.35 


8.8 


23.57 


4.55 


5.66 


0.45 


0.27 




2,4,7 


HD 18143C 


13.86 


8.66 


43.71 


5.2 


6.84 


0.28 


0.17 




2,4,7 


55 Cnc B 


13.15 


7.67 


79.8 


5.48 


7.21 


0.31 


0.20 




2,4,7 


HD 190360B 


14.4 


8.71 


62.92 


5.69 


7.71 


0.21 


0.10 




2,4,7 


Proxima Cen 


11.11 


4.38 


772.33 


6.73 


8.81 


0.21 


0.14 




2,4,7 


HIP 79431 


11.34 


6.589 


69.46 


4.751 


5.80 




0.35 


Jupiter 


2,8,9 


GJ 876 


10.16 


5.01 


213.28 


5.15 


6.65 




0.23 


2 Jupiters + Super-Earth 


2,8,10 


GJ 317 


12 


7.028 


109 


4.972 


7.22 




-0.20 


Jupiter 


2,11 


GJ 849 


10.42 


5.594 


109.94 


4.826 


5.80 




0.41 


Jupiter 


2,8,12 


GJ 179 


11.96 


6.942 


81.38 


5.018 


6.49 




0.20 


Jupiter 


2,8,13 


GJ 832 


8.66 


4.501 


201.87 


4.159 


6.03 




-0.24 


Jupiter 


2,8,14 


GJ 649 


9.7165 


5.624 


96.67 


4.0925 


5.55 




-0.03 


Saturn 


2,8,15 


GJ 436 


10.67 


6.073 


98.61 


4.597 


6.04 




0.10 


Neptune 


2,8,16 


GJ 581 


10.57 


5.836 


160.91 


4.734 


6.87 




-0.22 


Neptune + 3 Super-Earths 


2,8,17 


GJ 674 


9.36 


4.855 


220.24 


4.505 


6.57 




-0.24 


Neptune 


2,8,18 


GJ 176 


9.97 


5.607 


107.83 


4.363 


5.77 




0.06 


Neptune 


2,8,19 


GJ 1214 


14.71 


8.782 


77.2 


5.928 


8.22 




0.28 


Super-Earth 


2,20,21 
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